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Abstract. The nearby BL Lacertae (BL Lac) object Markarian 421 (Mkn 421) at a red shift z = 0.031 was observed 
to undergo strong TeV 7-ray outbursts in the observational periods from December 1999 until May 2001. The 
time averaged flux level F{E > 1 TeV) in the 1999/2000 season was (1.43 ± 0.04) TO -11 ph cm~ 2 s~\ whereas in 
the 2000/2001 season the average integral flux increased to (4.19±0.04) TO -11 ph cm~ 2 s _1 . Both energy spectra 
are curved and well fit by a power law with an exponential cut-off energy at 3.6(+0.4 — 0.3) a tat (+0.9 — 0.8) a y S TeV. 
The respective energy spectra averaged over each of the two time periods indicate a spectral hardening for the 
2000/2001 spectrum. The photon index changes from 2.39 ± 0.09 sta t for 1999/2000 to 2.19±0.02 sta t in 2000/2001. 
The energy spectra derived for different average flux levels ranging from 0.5 to 10 TO -11 ph cm s follow a 
clear correlation of photon index and flux level. Generally, the energy spectra are harder for high flux levels. From 
January to April 2001 Mkn 421 showed rapid variability (doubling time as short as 20 minutes), accompanied with 
a spectral hardening with increasing flux level within individual nights. For two successive nights (MJD 51989- 
51991, March 21-23,2001), this correlation of spectral hardness and change in flux has been observed within a 
few hours. The cut-off energy for the Mkn 421 TeV spectrum remains within the errors constant for the different 
flux levels and differs by AE = 2.6 ± 0.6 s tat ± 0.6 sys TeV from the value determined for Mkn 501. This indicates 
that the observed exponential cut-off in the energy spectrum of Mkn 421 is not solely caused by absorption of 
mufti- TeV photons by pair-production processes with photons of the extragalactic near/mid infrared background 
radiation. 
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1. Introduction 

The nearby BL Lac object Mkn 421 (z = 0.031) was the 
first extra-galactic source of very-high e nergy gamma-ray 
emission detected ( Punch et al., 199^ ). The source was 
monitored closely with the HEGRA imaging Cherenkov 
telescopes ( Konopelko et al., 1999 ) during the observa- 
tional periods from 1997-1998 while the source remained 
in a low flux state with an average flux level of a third 
of the flux observed from the Crab-Nebula and a power 
law energy spectrum with a photon index of 3.09 ± 0.07 
( Aharonian et al., 1999b| ) . It is well known that the source 
exhibits time variability on a sub-hour timescale at TeV- 



Table 1. Individual observation times, and mean zenith 
angles (9) are listed for the different observational periods. 



energies (Gaidos et al., 1996), which constrains the size of 
the emission region (R < 10 15 cm(5/10) (T var /lhr)) and 
indicates a relativistic bulk motion of the emitting plasma 
with a Doppler factor S = (T - ^T 2 - lcostf)- 1 > 10 at 
an angle 9 to the line of sight to avoid severe internal 



absorption (Celotti et al., 1998 

Th e observed variations of the source in the optical 



marginally with the TeV flux ( Nied et al., 2001 ; Jordan 
et al., 2001 ). However, the X-ray activity shows good cor- 
relation with the TeV flux level. During the observational 
periods December 1999 until May 2001, extended multi- 
wavelength observation campaigns were conducted by X- 
ray satellites (RXTE, BeppoSAX) and TeV observatories 
(HEGRA, CAT and VERITAS collaborations). During co- 
ordinated simultaneous observations in January/February 
2001 with RXTE, Mkn 421 showed a soft X-ray spectrum 
consistent with a peak position of the synchrotron spec- 
trum below 3 keV. The X-ray flux integrated between 
2-10 keV reached a value of 1.6 ■ 10 -9 erg/(cm 2 s) on 



January 31, 2001 ( [Horns et al., 200 1| ) where Mkn 421 
showed also a high TeV flux level. This is the highest X-ray 
state for Mkn 421 reported during a pointed observation. 
Interpretation of the broadband spectral energy distribu- 
tion (SED) in the framework of synchrotron-self-Compton 
models is discussed in Sect. [|. 

The TeV energy spectrum of Mkn 421 has been mea- 
sured by several instruments (HE GRA, Whipple, CAT) at 
different flux levels of the source ( Aharonian et al., 1999b ; 
Krennrich et al., 1999| , |200l| ; |Piron et al., 200l|) . The~ob- 
served spectral shape shows tentative indications for vari- 
ations and a possible correlation with the overall flux level. 
However, since the instruments have largely different en- 
ergy thresholds and systematic uncertainties, a firm claim 
of spectral variations of the TeV emission from Mkn 421 
could so far not be made. For Mkn 501, observations of a 
spectral hardening during the 1997 outbursts in TeV emis- 
sion have been claimed ( Djannati-Atai et al., 199S), but 
could not be verified by other experiments (Aharonian et 
al., 2001b| ). The TeV-spectrum of Mkn 501 softened during 
the subsequent years following the strong outburst while 
the mean flux has decreased by a factor of 10 to 1/3 of 



Date 


Year 


Obs.Time 


(9) 






\hrs~\ 


L J 


December 3-11 


1999 


4.8 


17.8 


January 31-February 15 


2000 


58.4 


20.9 




2000 


17.9 


13.2 


March 24-March 29 


2000 


28.4 


19.8 


April 23-May 5 


2000 


21.6 


17.1 


November 24-December 6 


2000 


20.8 


30.5 


December 28 


2000 


2.0 


16.4 


January 17- February 5 


2001 


94.2 


21.2 


February 13-February 28 


2001 


63.2 


20.5 


March 17-March 31 


2001 


51.6 


17.8 


April 14- April 26 


2001 


8.1 


20.6 


May 9- May 23 


2001 


14.9 


20.7 


Total 


385.9 


20.3 



the flux observed from the Crab-Nebula ( Aharonian et al 



and ra lio energy region have been found to correlate only 2001b ) 



Mkn 42 1 exhibited strong activity from January to 

I) with a peak diurnal flux 



May 2001 flBorst et al., 2001 



level on March 23/24, 2001 (MJD 51991.9-51992.2) with 
a value of (12.5±0.4) -lO" 11 ph cm~ 2 s" 1 above 1 TeV cor- 
responding to 7.4 times the flux of the Crab-Nebula. Only 
6 days before reaching the peak flux, the diurnal flux had 
been as low as (0.5 ± 0.1) -10 -11 ph cm~ 2 s" 1 , a factor of 
25 lower. 

We report on observations of Mkn 421 carried out with 
the H EGR A syst em of 5 imaging air Cherenkov telescopes. 
Sects. 3.1 & 3.2 contain the results on flux and spectral 
shape of Mkn 421 during the observational periods from 
December 1999 until May 2001. We show clear evidence for 
variations of the spectral shape of a TeV 7-ray source, get- 
ting harder with increased flux level (Sect. |3.3| ). Intranight 
flux and spectral variations are evident in at least two 
nights of observations (Sects. 3.4 & 3.5). The difference 
between cut-off energies present in the TeV spectra of 
Mkn 421 and Mkn 501 is discussed in Sect. 
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2. Data selection and analyses 

The data used here were collected from December 1999 
until May 2001 and comprises in total 385.9 hours of ob- 
servations under stable weather conditions, good atmo- 
spheric transparency and good detector performance (see 
Table [l] for a breakdown of the observation time for indi- 
vidual months). Data taken up to zenith angles of 45° are 
used which constitute more than 95 % of the data taken 
in total. 

Based upon comparing the observed rate with the ex- 
pected cosmic ray event rate, individual runs were rejected 
if the relative deviation from the expectation exceeds 20 %. 
In total 35 hrs of data fail to pass this selection. 
A large fraction of data were taken with a 5-telescope 
setup. However, in January 2000 and April 2001, only 
4 telescopes were operating. The collection areas for the 
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Fig. 1. For the observational periods 1999/2000 and 2000/2001 the diurnal integral fluxes observed above 1 TeV are 
indicated. The highlighted days (circled) are subject to separate discussion in Sect. 3.4fc3.5. The bottom part of each 
panel shows the hardness ratio calculated by integrating the flux in two separate energy bands covering 0.75 -1.5 TeV 
and 1.5 -4.0 TeV respectively. For short diurnal exposures, consecutive measurements have been combined. 
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Table 2. Listed arc the observational time (T & s ), the integral flux above 1 TeV and the results of x 2 minimisation 
on the differential photon energy spectrum using a function of the form dN/dE = N (E/TeV)~ a exp(— E/E ). The 
cut-off energy E has been fixed to the value derived from a fit to the combined data set taken from December 2000 
until May 2001. The quoted errors are statistical errors only. The lower part of the Table summarises the result of 
applying fits of a power law function to the time averaged spectra and examples of diurnal energy spectra with spectral 
variations during flares. 



Period 


T bs 


Fe>1 ToV 


N 


Q 


Eo 


X? ed (d.o./.)0 




[hrs] 




"lO" 11 " 




lO" 11 




[TeV] 








cm 2 s 




cm 2 s TcV 






dN/dE 


= N E- a exp(-£/3.6TeV) 












Fixed _Eq 





Dec 1999-May 2000 131.1 1.43 ±0.05 4.3 ± 0.3 2.39 ± 0.09 3.6 0.8(8) 

Nov 2000-May 2001 254.8 4.19 ±0.04 11.4 ±0.3 2.19 ± 0.02 3.6 1.2(8) 

0<F_ii<lf] 42.9 0.53 ±0.04 2.2 ±0.3 3.0 ± 0.2 3.6 0.3(3) 

KF_n<2 77.4 1.55 ±0.04 5.2 ±0.2 2.47 ± 0.08 3.6 0.9(6) 

2<F-n< 4 54.1 2.76 ±0.05 8.8 ± 0.3 2.36 ± 0.05 3.6 0.9(6) 

4<F_n< 8 129.4 4.22 ±0.03 12.1 ±0.2 2.18 ± 0.02 3.6 0.8(6) 

8<F_n<16 17.8 10.3 ±0.1 27.3 ±0.5 2.06 ± 0.03 3.6 0.5(6) 

MJD51964/51965 3.6 10.6 ± 0.4 32 ± 1 2.36 ± 0.06 3.6 0.9(6) 

MJD51991/51992 4.1 12.5 ± 0.4 31 ± 1 2.04 ± 0.05 3.6 1.2(6) 



dN/dE = N E~ 



no cut-of 



Dec 1999-May 2000 






3.8 ±0.2 


3.19 ± 0.04 


2.9(8) 


Nov 2000-May 2001 






10.1 ±0.1 


3.05 ±0.01 


53(8) 


MJD51989/51990,preflare 


1.9 


3.6 ±0.3 


8.1 ±0.7 


2.9 ± 0.1 


0.5(6) 


MJD51989/51990,flare 


1.1 


7.54 ±0.05 


12.0 ± 1.0 


2.2 ± 0.1 


1.0(5) 


MJD51990/51991,preflare 


1.5 


3.1 ±0.4 


8.3 ±0.8 


3.2 ± 0.2 


0.4(5) 


MJD51990/51991,flare 


1.0 


10.0 ±0.7 


15.2 ±0.9 


2.5 ±0.1 


0.7(5) 



Xred 

F-u 



X 2 /d-o.f. 

F(E > lTeV)/10" 



L phcm 



spectral analysis of data taken with a 4 or 5 telescope setup 
are calculated separately to take the difference in accep- 
tance between the two setups into account. Additionally, 
the degrading mirror reflectivity and the conversion fac- 
tors relating the registered digitised Cherenkov amplitude 
values in the cameras to the number of photoelectrons 
detected have been determined for monthly observational 
periods and are applied to the Monte-Carlo simulation to 
calculate the collection areas for each period separately. 

All observations were carried out in the so-called wob- 
ble-mode. During wobble-mode observations the source di- 
rection is positioned ±0?5 shifted in declination with re- 
spect to the centre of the field of view of the camera. This 
observation mode allows for simultaneous estimate of the 
background (OFF) rate induced by charged cosmic rays. 
The analysis uses an extended OFF-region to reduce the 
statistical error on the background estimate. A ring seg- 
ment (180° opening angle), from 0?3 to 0?7 distance from 
the camera centre at the opposite side of the ON region 
has been chosen. The events suitable for spectral analysis 



are selected in a similar way as described in Aharonian et 



al. (19 )9b), applying only loose cuts on the image shape 
and the shower direction in order to minimise systematic 
effects from energy dependent cut efficiencies. The loose 
directional and shape cuts accept events with the recon- 
structed direction being within an angular distance of 0?22 



to the source and the image shape parameter mean scaled 
width mscw < 1.2. 80% of the registered 7-events pass 
both the directional and the shape cut. To ensure that 
the images are not truncated because of the limited field 
of view, images with a distance > 1?7 to the centre of the 
camera are rejected. At the same time the reconstructed 
core position is required to be within 200 m distance to the 
central telescope. The energy reconstruction requires two 
images with more than 40 photoelectrons after applying 
a two-stage tail-cut removing night sky background con- 
tamination of the image. The angle subtended between 
the major axes of at least two of the registered images 
is required to be larger than 20°. The number of excess 
events after applying all selection criteria used in this anal- 
ysis amounts to ~ 40 000, constituting the largest photon 
sample collected from a single source with the HEGRA 
telescopes. 

2.1. Spectral analysis method 

The energy reconstruction and the reconstruction of en- 
ergy spectra have been described in detail elsewhere 
( Aharonian et al., 1999a ). For a large part of the data, con- 
temporaneous observations of the Crab-Nebula^ are avail- 



1 The Crab-Nebula is frequently used as a standard candle 
for TeV astronomy 
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Fig. 2. Separated into the two observational periods 
(1999/2000 and 2000/2001), the time averaged spectra 
show differences in the slope of the spectrum being harder 
in 2000/2001 coinciding with an increase of the flux by 
a factor of 2.9 ± 0.1. The shaded areas indicate the sys- 
tematic errors on the spectral measurement, being most 
prominent in the threshold region. Both spectra are well 
fit by a power law with an exponential cut-off energy at 
3.6 TeV. The photon index changes by Aa = 0.17 ± 0.07 
as indicated by a fit to the ratio of the two energy spectra 
which is well described by a pure power law. The dashed 
line indicates a fit of a function as given in Eqn. | with 
c = (0.04 ± 0.07) TeV -1 indicating that E is consistent 
between the two measurements. 



able, allowing to compare and to verify the expectation 
for detection rates and cut-efficiencies for 7-induced air 
showers derived from Monte-Carlo simulations with data. 

In the possible presence of a deviation from a pure 
power law in the form of an exponential cut-off, the re- 
construction of the spectral shape beyond this cut-off is 
influenced by the energy resolution. Even with an event- 
by-event energy reconstruction accuracy of AE/E sa 20 % 
as achieved with the stereoscopic reconstruction technique 
of air showers, a considerable fraction of events with over- 
estimated energies contribute to the flux of neighbouring 
bins. 

This becomes increasingly important for energies well 
in excess of the cut-off energy. Different methods like de- 
convolution using suitable algorithms and forward-folding 
techniques are applicable and have been pursued. An ad- 
vanced energy reconstruction technique with an improved 
energy resolution of « 10% (Hofmann et al., 200C) has 
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Fig. 3. For three distinct flux intervals (H,M,L, see 
Table |^) , the energy spectra have been extracted and the 
flux ratios calculated. The energy spectra get considerably 
softer with decreasing average flux. The cut-off energy re- 
mains unchanged within the statistical error. 



been proven to show similar re sults as obtained with th e 
analysis technique used here ( Aharonian et al., 2001a ). 
However, the results presented here are obtained with the 
conventional energy reconstruction method being less sen- 
sitive to changes in the detector performance and therefore 
less susceptible to systematic errors. 

The influence of the energy resolution has been ab- 
sorbed in the collection area A(E leco , 9), 9 indicating the 
zenith angle. The method applied here follows the same 
approach as described in Aharonian et at. (1999a), us- 
ing a collection area A 



A(E leco ,,6) being a function of 
the reconstructed energy E Teco . The collection area there- 
fore depends upon the spectral shape assumed. An it- 
erative method is applied to test a spectral hypothesis 
on the data. The procedure starts with an assumed arbi- 
trary spectral shape. For the sake of simplicity, a Crab-like 
power law spectrum with a photon index of 2.6 is assumed. 
If the reconstructed spectrum shows deviations from this 
shape, a new hypothesis can be tested by using a different 
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Fig. 4. The energy spectra get harder with increasing flux 
as is clearly seen when splitting data in five separate flux 
intervals and fitting a power law with a fixed exponential 
cut-off energy to the individual energy spectra. The indi- 
cated errors on the photon index are the statistical errors. 
The solid and dashed lines indicate fits of polynomial func- 
tions of the logarithm of the flux to the data (see text for 
details). Also indicated are the values for the night with 
the highest flux and as an open symbol the ar chival data 
from t he 199 7/1998 observations of Mkn 421 ( Aharonian 
et al., jL999b|) . 



spectral shape to derive the collection area and following 
the same procedure once more. The process is iterative 
and usually converges after two iterations. 

The method has been carefully checked with simulated 
energy spectra and also against other methods and has 
proven to be less dependent upon the exact modelling 
of the detector response with Monte-Carlo methods than 
more sophisticated procedures as mentioned above. 

2.2. Systematic uncertainties 

The systematic errors include as most important contri- 
butions the uncertainties in the conversion factors used 
to calculate the number of photo-electrons from the digi- 
tised pulses of the photo multiplier tubes, non-linearities 
of the electronic chain and the conservatively estimated 
variations of the response of the telescopes in the thresh- 
old region. Additionally, the uncertainties in the spectral 
shape enter the calculation of the collection area and cause 
an additional systematic uncertainty especially at energies 
well above the exponential cut-off energy. 

Besides the mentioned uncertainties, the overall flux 
determination is subject to a 15 % uncertainty in the cali- 
bration of the absolute energy scale of the telescopes. This 
uncertainty only affects the absolute flux calibration and 
does not influence the reconstructed shape of the spec- 
trum. The stability of the calibration of the energy scale 
has been carefully tested with data taken on the Crab- 
Nebula constraining possible time dependent effects seen 
as variations on the reconstructed flux to be less than 10% 



which translates into less than 6% uncertainty on changes 
of the energy scale (assuming a power law energy spectrum 
with a photon index of 2.6). A more in-depth discussion of 
the systematic errors is given in Aharonian et al. (1999a). 



3. Results 
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Fig. 5. a) The light curve includes the rise and decay 
of a flare observed during the night March 21/22, 2002. 
The solid line indicates a fit of a function F(t) as given in 
Eqn. The hardness ratio shows a hardening of the spec- 
trum correlated with the increased flux, b) The lower 
figure shows the fastest rise time observed in the data. A 
fit of an exponential profile (Eqn. ||) to the light curve 
indicates a doubling time of 21 ± 2 minutes and a decay 
time of ~ 4.5 hrs. However, the sensitivity to a specific 
profile of the flare is limited by the accuracy of the flux 
measurement. A Gaussian function (dashed curve) is also 
an acceptable description of the data. Again, as in Fig. ||a, 
the energy spectrum shows a hardening and softening cor- 
related with the increase and decrease of the flux. 
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3.1. Light curve 

The observed integral flux above an energy of 1 TcV 
for individual nights is shown in Fig. [l]. In the year 
2000, the highest diurnal flux level reached (6.4 ± 
0.3) TO" 11 ph cm _2 s _1 ( MJD 51662/51663). Note the 
smooth increase and decrease of the flux level in the course 
of a few days observed from MJD 51659 until MJD 51665 
(May 2000). The flux averaged from December 1999 to 
May 2000 is (1.43 ± 0.04 stat ± 0.2 sys ) -lO^ 11 ph cm" 2 s" 1 . 
In the observational period from November 2000 until 
May 2001 the average flux and the amplitude of vari- 
ations increased. The maximum diurnally averaged flux 
was observed in MJD 51991/51992 (March 23/24 2002) 
with (12.5 ± 0.4) TO- 11 ph cm" 2 s^see also Sect. p| ). 
The time averaged flux for the 2000/2001 measurements 
is (4.19±0.04 stat ±0.4 sys ) TO" 11 ph cm^s" 1 , which cor- 



E Q = 3.6 



responds to 2.4 times the flux of the Crab-Nebula (the 
integral flux above 1 TeV observed from the Crab-Nebula 
is (1.76±0.06 sta t±0.51 

iooop). 



•10 



-11 



ph cm 2 s 1 ( Aharonian 



et al. 



3.2. Time averaged energy spectra 

The observational periods 1999/2000 and 2000/2001 are 
treated separately for two reasons: The average flux level 
and the amplitude of variability is largely different for the 
two seasons (see Fig. [j]). Secondly, the two measurements 
are separated by one year and allow to investigate secu- 
lar changes by comparing the energy spectra for the two 
observational seasons 

A fit of a power law function dN/dE = N Q (E/ TcV)- a 
to the energy spectrum derived from the observations 
in the years 1999/2000 results in N = (3.8 ± 0.2 sta t ± 



0.4 S 



•10- 



ph cm" 2 s" 1 TeV" 



3.19 ± 0.04 stat ± 



0.04 sys with Xred(^-°-/-) = 2.9(8). This makes a power law 
fit seem unlikely given a probability of 5 • 10 -3 for a larger 
X 2 -value. The soft spectral shape is consistent with pre- 
vious measurements (Aharonian et al., 1999b). A fit of a 



power law with an exponential cut-off 



dN »r ( E 

= N • 

dE u V TeV 



E 



(1) 



results in a fit to the data with N = 4.3 ± 0.3 stat ± 
0.5 sys TO- 11 
0.04 sys , E = 



ph cm s~ 
3.8 (+*) 



0.9(7). The large statistical error on the position of the 
cut-off is a result of the strong correlation of the two pa- 
rameters (a,E ) used to characterise the spectral shape. 

The energy spectrum derived from observations in the 
period from November 2000 until May 2001 with a sub- 
stantially larger number of photons and therefore reduced 
statistical fluctuation is not consistent with a pure power 
law function (\ 2 — 424 with 8 degrees of freedom). A fit 
of a power law function with an exponential cut-off results 
in a smaller % 2 -value of 8.8 with 7 degrees of freedom and 
A^o = 11.4 ± 0.3 stat ± 0.4 sys , a = 2.19 ± 0.02 stat ± 0.04 sys , 



TeV. Table | shows the indi- 



h0.4\ (+0.9\ 
-0.3/gtat V-0.8/ S ys 

vidual results on the power law fit and fitting a power law 
with an exponential cut-off to different time periods and 
flux levels. 

The energy spectra derived individually from the two 
observational periods are presented in Fig. ||. The differ- 
ential flux is multiplied by E 2 to emphasise the subtle 
differences between the two spectraPJ Provided that the 
exponential cut-off is at the same energy for both years 
(Eq = 3.6 TeV), the ratio of the two spectra should follow 
a power law. The lower panel of Fig. |^ shows the ratio of 
the two energy spectra. Systematic effects that apply to 
both data sets in a similar manner cancel out for the ratio 
of the two spectra. As indicated with a solid line in the 
lower panel of Fig. ||, a power law fit is a good description 
of the ratio. The photon indices differ by Aa = 0.17±0.07 
indicating a trend for the spectr um to become harder at 
a higher flux level (see also Sect. . The dashed curve 
is the result of a fit with an exponential as expected for 
different cut-off energies: 

dN,, 



^ ' dE 



Ni 



h/h 

Aa 



E 



-A 



E 
TeV 
1 • exp( 
a 2 



■ exp(-E/Ei) 
-c-E) 



(2) 



c = E 1 

The result of a x 2_n t to /i / / 2 (see Fig. ||) results in a differ- 
ence c = (0.04±0.07) TeV- 1 (AE = (0.5±0.8) TeV) which 
is consistent with the two cut-off energies being identical. 

3.3. Correlation of flux and photon index 

To investigate possible correlations of the flux level and 
the spectral shape, the data have been split in five flux 
intervals as indicated in Table |^. The energy spectra for 
three of the flux bins (labelled Low, Medium, and High) 
are shown in Fig. |3|. The different flux levels show different 
spectral slopes. The ratio of the energy spectra is again 
well-described by a power law indicating that the cut-off 
energy for the different flux levels remains constant within 
the statistical error, whereas the photon index changes by 
Aa = 0.31 ± 0.04 for the ratio of the H and M. For an 
average flux of F-u < 1 (Low), the spectrum softens such 
that the spectral index changes by Aa = 0.9 ± 0.2 with 
respect to H. 

The general trend of hardening of the energy spectrum 
with increased flux is shown in Fig. || where the photon 
index a is shown as a function of integral flux above 1 TeV. 
To paramctcrise the dependence of photon index and flux, 
a function of the form 

2 

Ct(F E >l TcV) = X! ai l °Sw( F E>l TcV) 
i=0 



This corresponds to vF(y): 1 TeV/(cm 2 s)=1.6 erg/(cm 2 
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stant background is fit to the light curves: 



23:00 



0:00 



1:00 



2:00 



3:00 



UTC [h] 



Fig. 6. The lightcurve measured over the course of 
3.6 hrs shows the decaying part of a flare that proba- 
bly started before observations began. This light curve 
is well described by a simple exponential function 
F(t) = fa cxp(-t/T decay ) with a decaying time T decay = 
4.5(+0.5 — 0.4) hr. The hardness ratio (lower part of the 
figure) for this night remains constant within the errors 
and the low value of the hardness ratio indicates a soft 
spectrum despite the high flux level. 



a = -2.7±0.1 
ai = 1.1 ±0.3 
a 2 = -0.4 ±0.2 



with Xred(d-°-f ) = 1-1(2) represents the data points quite 
well. A simpler functional dependence of the form 

a(F E> i Tcv) = b Q + bi ■ log 10 (FE>i Tc v) 
b = 2.7 ±0.1 
&i = 0.47 ±0.02 

is compatible with the data with a x 2 ed (d.o./.) = 2.3(3) 
(indicated by the dashed curve in Fig. |I|). A % 2 -fit of a 
constant value is excluded by a very large X^ ed {d.o.f.) = 
14(4). 

It is interesting to note that the energ y spectrum of 



Mkn 4 21 observed during 1997 and 1998 ( Aharonian ct 
al., 19 j)9b| ) fits nicely in this picture when applying the 
same fitting method to the archival data keeping the cut- 
off energy fixed at 3.6 TeV (see Fig. ||). This data-point 
has been excluded from the fitting of the functional de- 
pendence of the photon index. 

3.4. Diurnal flux variations 

Intra-night variations of the flux are seen during a large 
fraction of the observational nights. In the past, doubling 



times as short as 15 minutes have been observed (Gaidos 



et al., 1996). To quantify the rise and decay times of flares, 
a simple profile consisting of two exponentials with a con- 



F(t) = f a + fi 



exp(T risc • (t - t )) t < t ,„s 

exp(-T dcC ay ' {t - t )) t > t 

An example of a night with a fast decay time (March 21/22 
2001, MJD 51989.93-51990.14) is shown in Fig. |a. The 
exponential decay time as given by the result of a x 2 -fit 
of a function F(t) (Eqn. ||) to the light-curve results in 
Tdecay = 15(+9 — 3) min, whereas the rise time turns out 
to be T r i SO = 46(+60— 15) min. This result is indicating an 
asymmetry of the light curve during a flare. However, the 
presence of a possible sub-structure in the light curve (at 
«1:00 UTC) which is not resolved could be responsible for 
an overestimated rise time. The night with the fastest rise 
time (March 22/23 2001, MJD 51990.925-51991.12) of the 
data set shows a smooth increase of the flux by a factor 
of 5 within 3 hours (see Fig. ||b). Fitting again a function 
F(t) to the light-curve results in T r i se = 0.52 ± 0.03 hr. 
This translates into a doubling time of 21 ± 2 minutes. 
The decaying part of the light curve shows indications for 
a decay time of 4.5 hr. However, a fit of a Gaussian func- 
tion with a constant background (dashed line in Fig ||b) 
describes the data equally well and would increase the es- 
timated doubling time to « 30 min. 

Another example of a diurnal light curve is shown 
in Fig. | (February 24/25 2001). The light curve of this 
night is well described by a simple exponential decay with 
Tdecay = 4.46(+0.54 — 0.39) hrs similar to the decay time 
seen in the night March 23/24 2001 (see Fig. |b). It is 
conceivable that the observed light curve in this particu- 
lar night covers the end of a flare for which the first part 
was not visible for HEGRA. Another peculiarity of this 
night is the small and constant hardness ratio indicating 
a soft spectrum at a comparably high flux level (see also 
next Sect.). 

3.5. Diurnal spectral variations 

Even with the large number of photons detected, the sen- 
sitivity for intra-night spectral changes is limited. In a few 
cases, the observational time and the number of detected 
7-rays is sufficient to perform spectroscopy on time scales 
of hours. A particular set of flares already mentioned in 
the previous Sect, occurred during the night MJD 51989.9- 
51990.2 (March 21/22 2001) (see Fig. |a) and the sub- 
sequent night MJD 51990.9-51991.2 (March 22/23 2001) 
(see Fig. ||b) . The light curves show in both cases a rising 
flux with a doubling time of less than 1 hr. At the end of 
the flare, the flux smoothly decreased. During the rising 
time of the flux, the hardness ratios exhibit in both cases 
a spectral hardening, well correlated with the increase of 
flux and a spectral softening at the end of the flare. For 
the night March 21/22 the hardness ratio as well as the 
flux level return after the flare to the initial state it had 
prior to the flare. For the night March 22/23, the spec- 
trum softens to the level it had before the flare, whereas 
the flux remains at a high level at a factor of 3 higher than 
before the flare. 
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Fig. 7. a)-b) The two nights from March 21/22 (MJD51989.9-51990.2) and March 22/23 (MJD51990.9-51991.12) 
show a spectral hardening while the flux increases. For both nights, preflare and flare spectra have been extracted and 
compared directly. In both cases, the photon index changes by Aa w 0.75 while the flux increases by a factor of 3-4. 



To gain further information on the spectral shape for 
the individual nights with indications for spectral varia- 
tions, energy spectra were extracted for two time inter- 
vals (indicated by the vertical lines in Fig Sa-b) from each 
night. The time intervals were chosen to cover the preflare 
and the flare state. For both flares, the preflare energy 
spectrum is well described by a power law with a photon 
index of rj 3.0, and a hardening by Aa = 0.75 ± 0.25 
during the flare (see Fig. Ra-b) . 

However, light curves with flux variations have been 
observed that are not accompanied by spectral variations. 
An example of a night with a varying flux and a constant 
spectral shape is given in Fig. ||a, where the energy spectra 
for two time intervals from the night MJD 51964/51965 
have been extracted (marked by vertical lines in Fig. ^). 
The spectrum is despite the large flux (F(E > 1 TeV) = 
(10.6 ± 0.4) -10~ n ph cm _2 s _1 ) comparably soft with a 
photon index of 2.85 ± 0.09 for the first half of the night 
and 2.95 ± 0.09 for the second half of the night. 

It is also interesting to note after comparing individ- 
ual nights, that the spectral shape does not always re- 
late directly to the absolute flux-level. As an example, 
Fig. ||b displays the spectrum of the night with the high- 
est diurnal flux (MJD 51991/51992) with a level only 20 % 
higher than the night February 24/25 discussed earlier. 
During this night, no strong indication for flux variability 
or changes in the hardness ratio are evident. The energy 
spectrum is hard with a photon index of a = 2.04 ± 0.05 



(Fig. ||b). The ratio of the two spectra indicate a difference 
of Aa = 0.36 ±0.09. 



4. Mkn 501 and Mkn 421 cut-off energies 

In the presence of extragalactic background light (EBL) at 
wavelengths of A = 0.5 ... 30 /im, TeV-energy photons are 
subject to absorption by pair-production processes with 
the l o w energy backgrou nd photons ( Gould fc Schreder 
1966 ; Stecker et al., 1992 ). As a result, the observed spec- 
trum of extragalactic TeV sources is modified. A very 



likely feature in the observed spectrum is an exponential 
cut-off. However, the presence of a cut-off in the observed 
spectra could also be caused by processes related to the 
production of TeV photons (e.g. decrease of the Compton- 
scattering cross-section for large centre-of-momentum en- 
ergies, absorption inherent to the source). In the simplest 
picture of a source spectrum that follows a power law, the 
position of the cut-off should be identical for all sources at 
the same red-shift provided that the EBL is isotropical. 

The two well-established extragalactic sources of TeV- 
radiation (Mkn 421 & Mkn 501) happen to almost coincide 
in red shift, which makes an investigation of differences in 
the position of the cut-off energies meaningful. 

In Fig. ^|the energy spectrum of Mkn 501 as measured 
during 1997 with a cut-off energy of 6.2 ± 0.4 stat (+2.9 - 
1.5) sys TeV and the energy spectrum derived from ob- 
servations of Mkn 421 in 2000/2001 are compared. The 
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Fig. 8. a) For the night February 24/25, two energy spectra at different flux level have been extracted. No indication 
for a spectral evolution is visible. Both energy spectra are well described by a pure power law with a soft photon 
index of 2.85 ± 0.09 and 2.95 ± 0.09 respectively, b) To compare two energy spectra at similar flux level but with 
different spectral shape, the energy spectra were extracted for the night with the highest diurnal flux (March 23/24, 
MJD 51991/51992) and compared with the diurnal spectrum from February 24/25. 



upper panel contains the differential energy spectra mul- 
tiplied by E 2 together with the respective fit functions 
indicated as solid lines. The dashed curve is the result of 
a fit of a power law with a fixed exponential cut-off at 
Eq — 6.2 TeV to the Mkn 421 data points. Ignoring possi- 
ble systematic errors, the resulting y^ eA {d.o.f.) = 5.7(12) 
excludes Eq = 6.2 TeV as a cut-off energy for the Mkn 421 
spectrum. The lower panel of Fig. || shows the ratio of 
the two energy spectra. A pure power law fit to the data 
is very unlikely because of a large ^ cA (d.o.f.) = 3.6(12) 
whereas a function as defined in Eqn. |^ results in a ac- 
ceptable Xr cd ( d -°-f-) = 1-5(12) with Aa = 0.35 ±0.08 and 
c = (0.12 ± 0.03) TeV -1 . We note that constraining the 
fit region to energies above 2 TeV where absorption starts 
to become important, a pure power law fit is not excluded 
any more for the ratio. 

There is an overlap of the cut-off energies determined 
for Mkn 421 (E =3.6(+0.4 - 0.3) sta t(+0.9 - 0.8) sys TeV) 
and Mkn 501 (Eg =6. 2 ± 0.4 stat (+2.9 - 1.5) sys TeV) 



(Aharonian et al., 1999a) taking the combined systematic 
and statistical errors into account. However, the origin of 
the systematic errors are very similar for both results and 
therefore strongly correlated. As mentioned in Sect 



2.2, a 



6 % systematic uncertainty on the relative energy scale re- 
mains when comparing the two energy spectra. Combining 
the independent statistical uncertainty and the estimated 



systematic uncertainty on the relative energy scale results 
in AE = 2.6 ± 0.6 sta t ± 0.6 sys TeV for the difference on 
the energy of the cut-off position. A detailed investigation 
of possible differences of energy spectra will be the issue 
of a forthcoming paper. 



5. Multi- wavelength observations 

The flux increase observed from Mkn 421 in the TeV en- 
ergy band has been accompanied by a strongly increased 
X-ray flux. Detailed studies of the correlation of the vari- 
ability in both energy bands are under-way. The simul- 
taneously taken data from January /Feb ruary 2001 with 
RXTE and HEGRA ( |Horns et al., 200l| ) reveal a largely 
different spectral evolution for Mkn 421 in the X-ray and 
TeV being almost opposite to what has been observed 
from Mkn 501: Whereas the shape of the TeV-spectrum 
of Mkn 501 remained constant during the 1997 activity 

the Mkn 421 spectrum 



period (Aharonian et al., 1999a), 
showed a spectral hardening with increasing flux. A dif- 
ferent spectral evolution is also seen in the X-ray spec- 
trum for the two sources: Whereas the X-ray spectrum of 
Mkn 501 is very hard during the flare which is commonly 
interpreted as a shift of the synchrotron peak position 



^ne.ak f rom below 1 keV to 100 keV and back (Tavecchio 



|st al., 2001 ), the X-ray spectrum of Mkn 421 remains soft 
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Fig. 9. Comparing the energy spectra of Mkn 501 
(Aharonian et al., 1999a) with the measurements of 



Mkn 421 (2000/2001): The observed differential spectra 
are well described by a power law with exponential cut- 
off (see Eqn. Both, the cut-off energy and the pho- 
ton index differ for the two energy spectra (oiMknbm — 
1.92±0.03 sta t±0.2 sys , a Affe „42i = 2.19±0.02 stat ±0.04 sys ). 
A fit of a fixed cut-off energy at 6.2 TeV as measured 
for Mkn 501 while letting the photon index vary to the 
Mkn 421 data results in the dashed curve in the upper 
panel with x 2 cd = 5.7(12), excluding the larger cut-off 
energy provided that systematic effects are negligible. To 
further study the difference of the energy spectra and to 
reduce the impact of possible systematic effects, the ratio 
of the two energy spectra are calculated. The flux ratio 
is not compatible with a pure power law model (dashed 
curve, Xrcd = 3.6(12)) indicating a difference in the shape 
of the two energy spectra. A function of the form given in 
Eqn H gives an acceptable fit (j£ eA [d.o.f.) — 1.5(12)) with 
Aa = 0.35 ± 0.07 and c = (0.12 ± 0.02) TeV- 1 . 



with a moderate shift of the synchrotron peak f rom be- 
low 1 keV at a flux level of ~ 10~ 10 erg/(cm 2 s) (Malizia 
et al., |2000D to w 2 - 3 keV on January 31,2001 when 
the flux measured between 2-10 keV reached -F2-10 = 
1.6 • 10 -9 erg/(cm 2 s) (Horns et al., 2001). These observa- 



tions confirm the energy dependence of the peak-position 



sy 

peak 



F° 



Sy 



claimed to follow a power law behaviour with E 
for Mkn 421. For Mkn 501 the relation of the peak position 
and flux follows E s ^ eak cx Fg y underlining the very differ- 



ent behaviour of the two extreme BL Lac objects (Fossati 



et al, >00C) 



The TeV data of Mkn 421 taken simultaneously with 
the X-ray observations in January/February 2001 indicate 
a hard spectrum with a photon index a w 2 that deviate 
substantially from the quiescent energy spectrum being 
much softer (a ss 3). The constraints on a one-zone SSC 
model of the source derived from these multi-wavelength 
observations are severe. The X-ray spectrum with a peak- 
energy < 3 keV constrains the electron-spectrum which 
is responsible for the synchrotron peak emission to have a 



low energy ("/peak = 1 .5 • 10 5 <5 50 ' 5 B^*). As a result, with- 



out invoking large S and a small magnetic field, the SSC 
model predicts a peak energy for the Compton-scattered 
spectrum below 1 TeV and the spectral shape is expected 
to be soft in the energy region from 1-20 TeV because of 
the Klein-Nishina suppression of the Compton-scattering 
cross-section. However, the observed TeV-spectra are very 
hard and a small magnetic field B < 0.22 G and a large 
Doppler-factor S > 50 are required to achieve an ade- 
quate fit of the TeV da ta with a SSC-model (see e.g. 
Krawczynski et al. (2001 )). The absorption of TeV pho- 



tons due to pair production with photons of the extra- 
galactic radiation field adds to the difficulties, because 
the source spectra corrected for the absorption are intrin- 
sically harder than the observed spectra which are not 
easily described with a simple SSC-one-zone model (see 
e.g. [Aharonian (200l|)). 



6. Summary and conclusions 

The observations of Mkn 421 at TeV-energies indicate sig- 
nificant spectral variations for different flux levels. The 
photon index a changes from « 3 to f=a 2 with a flux in- 
creasing from F(E > 1 TeV) = 0.5 • F_n ... 10 • F_u. 
Diurnal flux variations have been observed with rise and 
decay times as short as 20 min. During particular nights 
where a complete flare was observed, the photon index 
determined during the preflare and the flare period differ 
by as much as Aa = 0.75 ± 0.25 getting harder while 
the flux increases. All energy spectra are well-fit by a 
power law with an exponential cut-off at Eg =3.6(+0.4 — 
0.3) s tat(+0.9 — 0.8) sys TeV without indications for varia- 
tions of the cut-off energy at different flux-levels. However, 
the sensitivity on changes of the cut-off energy is lim- 
ited by strong correlations between the two parameters 
(a, Eh) used to describe the energy spectra. A compari- 
son of E for Mkn 501 and Mkn 421 indicates that the 
position of the cut-off energy for Mkn 501 is larger by 
AE = 2.6 ± 0.6 sta t ± 0.6 sys TeV taking the combined sys- 
tematic and statistical error into account. A difference of 
the cut-off energies for two sources at similar red shift 
precludes the possibility of the cut-off being only an ab- 
sorption feature. 
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